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We derive the relativistic quantum transport equation for the pion dis- 
tribution function based on an effective Lagrangian of the QHD-II model. 
The closed time-path Green's function technique, the semi-classical, quasi- 



, particle and Born approximation are employed in the derivation. Both the 

> . 

I/"") mean field and collision term are derived from the same Lagrangian and 

' presented analytically. The dynamical equation for the pions is consistent 

m ; 

with that for the nucleons and deltas which we developed before. Thus, we 



obtain a relativistic transport model which describes the hadronic matter 
with N, A and tt degrees of freedom simultaneously. Within this approach, 
we investigate the medium effects on the pion dispersion relation as well as 
the pion absorption and pion production channels in cold nuclear matter. 
In contrast to the results of the non-relativistic model, the pion dispersion 



X 

relation becomes harder at low momenta and softer at high momenta as 
compared to the free one, which is mainly caused by the relativistic kinet- 
ics. The theoretically predicted free ttN — ► A cross section is in agreement 
with the experimental data. Medium effects on the irN — > A cross section 
and momentum-dependent A-decay width are shown to be substantial. 

PACS number(s): 24.10.Cn; 13.75.Cs; 21.65. +f; 25.70.-z 
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I. INTRODUCTION 

Pion physics is an important topic in nuclear physics. Recently, it received renewed 
interests in relativistic heavy-ion collisions because pions are the most abundantly pro- 
duced particles at relativistic energies. Studies of pionic many-body degrees of freedom in 
high-energy nucleus-nucleus collisions were initiated by Gyulassy and Greiner PJ and by 
Migdal Q. Since then, considerable efforts from both experimental |§- 0] and theoretical 



T0|j - |26| groups were made to study various aspects of the in- medium pion dispersion re- 
lation and pion dynamics, such as pion spectrum, pion and anti-pion flow in hot and dense 
nuclear matter. Due to the high interaction cross section of the pion with the nuclear 
environment they are continuously absorbed by forming A-resonances which then decay 
again into pions. Therefore, pions have a chance to be emitted during the whole course of 
the reaction. While the high-energy tail of the pion spectrum provides information about 
compressed and excited nuclear matter in the early reaction stage, the low-energy part of 
the pion spectrum and pion flow contain information of the in-medium pion potential and 
nuclear equation of state (EOS) [pi] ) p6(| . The low- and high-energy pions originate from 



different stages of the collision. A detailed study of the pion dynamics allows to extract 
the time evolution of heavy-ion collisions. 

On the other hand, dileptons produced from 7r + — 7r~ annihilation [p!7|,^7|-|30| provide 
information on the high density phase at time scales of 1 fm/c. Since dileptons can leave 
the reaction volume essentially undistorted by final-state interactions, as was first pointed 
out by Gale and Kapusta ||27|| , they are expected to be a good tool for an investigation 



of the violent phases of high-energy heavy-ion collisions. Recent data by the CERES 



collaboration BTJ] show a substantial modification of the dileptons yield which might be 



explained either by many-body effects [32] or by an enhanced p-meson production (via 



7r + — 7r~ annihilation) and a dropping p-mass in the medium p3| , |3^| . Indeed, the property 
of the p-meson as well as the A resonance are strongly influenced by the change of pion 
property in the medium due to the large p — > 7r + 7r~ and A — > Ntt decay widths. A detailed 
knowledge of pion dynamics in heavy-ion collisions is a prerequisite for a quantitative 
description of dilepton production at SIS and SPS energies. 



It was recently proposed that the difference between 7r~- and 7r + -spectra can be at- 
tributed to the influence of isospin and Coulomb fields p5| . This should allow to extract 
the effective Coulomb field at the instant of the average pion emission. The comparison 
of spectra of positively and negatively charged pions can also be used to learn about the 
freeze-out of the pions during the expansion phase PE| , |3"T|| . It then provides a method to 
determine the size of fireball during the nuclear expansion process. 

Since the importance of pions in heavy-ion collisions has been recognized for more 
than two decades one may believe that the elementary pion properties in the hot and 
dense nuclear matter are already well understood. Unfortunately, the situation is quite 
different to this expectation: understanding the pion dynamics in high-energy nucleus- 
nucleus collisions is still a major challenge to modern nuclear physics. This paradoxical 
circumstance seems to be mainly due to theoretical rather than experimental inadequa- 
cies. Experiments have been able to record pion spectrum |6|,|8|,j38|1 and pion flow and 
antiflow P,|3"5| with rather high accuracy. A reasonably strict treatment of pions in the 
transport theories is, however, still not available. Most theoretical approaches included 
the interaction of the pions with the surrounding nuclear medium only by collision pro- 
cesses. A free-particle assumption was usually assigned to pions, while it is well known 
that the pion dispersion relation will be changed substantially in the medium due to the 
strong p-wave interaction. Some authors pl| , p2| , p6| implemented the real part of the pion 
optical potential from the non-relativistic A-hole model [^D[ to study the pion spectrum. 
In Ref. p6| a phenomenological parameterization suggested by Gale and Kapusta p7| 
was also tried. However, different model treatments gave rather different results. None of 
the currently available models are able to reproduce the experimental spectrum over the 
entire range of energy. The source of the problem seems to be that the observed quanti- 
ties are sensitive to several of the unknown pion properties in the hot and dense matter. 
The most important ones are the in-medium pion dispersion relation (the real part of the 
pion self-energy) and the in-medium pion cross sections (the imaginary part of the pion 
self-energy). A self-consistent treatment of both the real part and imaginary part of the 
pion self-energy is necessary to obtain useful information from experimental observable. 
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However, in transport theories one usually uses the experimentally determined free cross 
sections and incorporates a dispersion relation from the non-relativistic A-hole model 
or simply employs the free-particle assumption. A self-consistent description of pions in 
transport models has not even been reached in the non-relativistic case. In the present 
work, we still don't have a fully self-consistent treatment. But we go a step further than 
the present available relativistic transport theories, i.e., we derive the in-medium pion 
dispersion relation and the in-medium pion cross sections from the same Lagrangian and 
treat them in a relativistic description. 

Another daunting obstacle to a quantitative description of pions in heavy-ion collisions 
is the width and shape of the A. As has been pointed out before, when a pion in the hot 
and dense matter collides with a nucleon, it will be absorbed to create a A. Then the A 
decays again into a pion- nucleon " pair" . Thus the amplitude for creating and absorbing 
pions will be sensitive to the in-medium A-decay width, which must be modified by the 
presence of matter due to the potential energies of N, A and tt. However, in the presently 
available transport models, the free A-decay width is commonly employed. From the 
theoretical point of view, realistic models for describing pions in dynamical processes 
should in principle at least also treat deltas and nucleons simultaneously in an unified 
framework. 

It is the purpose of this paper to develop the relativistic transport theory for pions 
within the framework of the relativistic Boltzmann-Nordheim-Vlasov-Uehling-Uhlenbeck 
(RVUU/RBUU) equation. The RVUU model has been applied successfully in studying of 
high-energy heavy- ion collisions [15|,|l]-^5]. In Ref. |15| we briefly discussed the possible 
extension to include the pion degree of freedom. By means of the density matrix method 
Wang et al. |16| developed a transport theory for the N, A and n system. In their work 
pions are treated as a free particle, the detailed expressions of the collision term are 
not given explicitly. On the other hand, in Refs. |47|]- |kJ we developed a set of self- 
consistent equations for N, A and iV*(1440) distribution functions in which both mean 
field and collision term are derived from the same effective Lagrangian and expressed 
analytically. However, mesons (a, u, tt) were treated as virtual particles. In a physically 
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reasonable scenario, the creation and destruction of real as well as virtual mesons ought 
to be described simultaneously and on the same self-consistent footing [1531 . Hence, one is 



forced to solve coupled Boltzmann equations not only for baryons but also for all relevant 
mesons. This will cause significant numerical difficulties and might be beyond the ability 
of modern computers. As a first practical step, let us here treat the pions explicitly. The 
pion is the most frequently observed meson. The other mesons still remain treated as 
virtual mesons. This is the main strategy of our present work. Here we should note that 
the transport equations for a and us mesons were discussed in Ref. |15| , |55|j based on the 



Walecka model |5q| . In Ref. [|5l| the mesons turn out to be treated as free particles due 
to the approximations used in that work. Also, no concrete expressions for the collision 
term were given there. We will come back to this point in Sect. III. 

Starting from an effective Lagrangian of the QHD-II model 5(| we here derive a RVUU 



equation for the pion distribution function in which both the mean field and the collision 
term are derived simultaneously and expressed analytically. In our framework a fully 
relativistic treatment is realized and medium effects are included. Furthermore, we treat 
N, A, and 7r in an unified framework based on the same effective Lagrangian and finally 
obtain a set of coupled equations for hadronic matter. The paper is organized as follows: 
in Sect. II we briefly review the closed time-path Green's function technique which plays 
a central role in our derivation. An effective Lagrangian for the N, A and it system 
interacting through the exchange of virtual mesons is also presented there. In Sect. Ill 
we derive the RVUU-type transport equation for the pion distribution function. The main 
ingredients of the equation are the relativistic mean field and collision terms, which are 
calculated from the same effective Lagrangian and presented analytically in Sect. IV and 
V, respectively. In Sect. VI we present the numerical results for the in-medium pion 
dispersion relation and A-formation cross section. Finally, the summary and outlook is 
given in Sect. VII. 

II. PRELIMINARIES 

In the present work we employ the closed time-path Green's function technique. For 
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a detailed description of this Green's function technique for non-equilibrium system, we 
refer to Refs. Here we give a brief review for the reader's convenience. In the 

Heisenberg picture the Green's function Gp(l, 2) of fermions and A#(i, 2) of bosons can 
be defined on the time contour depicted in Fig. 1 as 

iG F {l,2) = {T[* H {l)* H (2)]), (1) 
zA B (l,2) = (T[$h(1)$ h (2)]> - <$ H (1))<$ H (2)>, (2) 

where 1, 2 denote x\, X2; ^#(1) and ^#(2) represent the field operators of the nucleon and 
delta in the Heisenberg picture and $#(1) and $#(2) are those of the a, lu, it and p. Here 
we have specified the initial state by assuming that its density operator commutes with 
the particle- number operator ||5711 . Furthermore, we assume that the initial state admits 
the Wick decomposition (is noncorrelated) . Thus, in Eq. (1) the expectation value of a 
single fermionic field vanishes. In the case of bosonic Green's functions, the contributions 
from classical expectation values have been subtracted in order to concentrate on the field 
fluctuations around the classical values. On the other hand, the second term on the right- 
hand side of Eq. (2) explicitly indicates the presence of the mean field. According to the 
position of field operators on the time contour, we have four different Green's functions 
for fermions 

iG$+(l,2) = {T a S> H (l)* H (2)), 
<GJ-(1 > 2) = (* h (1)*h(2)>, 

iG F +(l,2) = -(* H (2)* H (1)), (3) 



and four for bosons 



iAi~(l,2) = <r c $ H (l)$tf(2)> - (^(1))($h(2)>, 
•A+t(l,2) = (n H (l)^ H (2)} - <$tf(l)><$tf(2)>, 
zA+-(l,2) = (d> H (l)<M2)> - ($ h (1)><$h(2)>, 

*A B +(1,2) = (d>„(2)<Ml)) - ($ H (1)>($ H (2)>. (4) 
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Here T c is the chronological ordering operator, T a is the anti-chronological ordering op- 
erator. The designations — and + are attributed to the respective time path shown in 
Fig. 1. 



Fig. 1 

We further on express the Gf(1, 2) and A#(l, 2) in a compact matrix form 

w ,=ir 2, ? 2) i 

iGp (1,2) *G++(1,2) 



and 

iA B (l,2) 



(6) 



<A+-(1,2) iA+ + (l,2) y 

It should be pointed out that the four Green's functions in Eq. (5) are not independent. 
They satisfy the following relations 

iGp~(\, 2) = 9(t x - t 2 )iG%-(l, 2) + 9(t 2 - tJiGp+il, 2), (7) 
iG F + (l, 2) = e(h - t 2 )iG F + (l, 2) + 6(t 2 - h)iG+-{l, 2). (8) 

Here 9(t\ — t 2 ) is defined as 

!1 U is later on a contour than t 2 
(9) 
t\ is earlier on a contour than t 2 

The same relations hold for the boson Green's functions in Eq. (6). 

In order to use the powerful perturbation expansion method of field theory, we choose 
the interaction picture. The time-ordered products in Eqs. (1) and (2) can then be 
rewritten as 

(T[* H (1)V H (2)]) = (T[exp(-i I tfeiZj(a;))ttj(l)*j(2)]), (10) 
<T[$ h (1)$h(2)]> = (T[exp(-i f dxH^x))® /(l)$/(2)]), (11) 
($„(!)) = (T[exp(-i I )$/(!)]). (12) 



Here ipi(l), '0/(2) and $/(l), $/(2) represent the field operators in the interaction picture; 
/ dx =/ dtdx, / stands for an integral along the time axis given in Fig. 1. The definition 
of Eqs. (3) and (4) and the relations of Eqs. (7), (8) are still valid in the interaction picture 
for both the full Green's functions Gf(1,2), A#(l,2) and zeroth-order Green's functions 
(i.e., non-interacting Green's functions) G%(1,2), A° B (1,2). The detailed expressions of 
the zeroth-order Green's functions as well as Hi in Eqs. (10) - (12) are determined by the 
specific effective Lagrangian used in the model. As a preliminary step towards a complete 
description of hadronic matter, we first consider a system consisting of real nucleons, 
deltas and pions interacting through the exchange of virtual cr, uo, tc and p mesons. In 
order to avoid extensive cancelations of large terms to correctly describe the small S-wave 
ttN scattering length, we choose the phenomenological pseudovector form for the 71NN 
and 7rAA coupling. With this choice of coupling, the value of the S-wave ttN scattering 
length turns out to be —0.010 |5(J while the empirical value is —0.015 ± 0.015 [59[|. The 



inclusion of the p-meson degree of freedom is important for the tttt scattering due to vector 
meson dominance |)(J . We furthermore include two non-linear meson coupling terms aim 
and pirn which are applied only to the tctt scattering. The total effective Lagrangian can 
then be written as 

C = C F + C h (13) 
Here Lp is the Lagrangian density for free nucleon, delta, and meson fields 

C F = #7^ - m nW + - m aWa 

+ ^ad»a - U(a) - + U(u) 

+ l -{d^d^ - mW) - + l -m\ Pil ■ p» (14) 

and U(cr), U(uj) are the self-interaction part of the scalar field []6T| and vector field [p2| , |63 



U(a) = l -mla 2 + ^b(g% N af + \c(g° NN a)\ (15) 
U(uj) = \ m l^(l + ^fl^), (16) 
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respectively. Here the field tensor for the rho and omega are given in terms of their 
potential fields by 



(17) 



and 



(18) 



The interaction Lagrangian Cj consists of baryon-baryon, baryon-meson and meson-meson 
terms, which are given by 



+gNNi>(xh»l5T ■ 1p(x)d^7T(x) - ^g P NN ^(x)^T ■ 1p(x)p^(x) 

+g AA ipAu(x)ip A (x)a(x) - g" J AA i) Av (x)-i^ A (x)uJ ,1 (x) 
+g AA i>A„{x)l^T ■ ^(x)^Tr(x) - ^ A <M*)7„T ■ r^P^x) 
-g AN i>^{x)d»Tz{x) ■ S+^(x) - gl N ^(x)S^(x) ■ d»ir(x) 
+ ^g^m a (y(x)^(x) ■ ir(x) + g (m [d^Tz(x) x n(x)] ■ p^x) 
= g^ N ^{x)Y^{x)^ A {x)+gi A ^ Au {x)Y^{^A{x) 

-g AN i>^{x)d»Tz{x) ■ S + ^j(x) - g AN i>(x)Si> Alx (x) ■ d»*{x) 
+gt*i{x)T\<K j {x)<$> A {x) 



(19) 



In the above expressions ip(x) is the Dirac spinor of the nucleon and ipA^x) is the Rarita- 
Schwinger spinor of the A-baryon. r is the isospin operator of the nucleon and T is 
the isospin operator of the A. S and S + are the isospin transition operator between 



the isospin 1/2 and 3/2 fields. g]^ N = j^/m^, g 



> BAN 



f*/m w ; = 7 A r A , T% = >y A T A: 



T A = "y A T A , A— a, cu, n, p, the symbols and notation are given in Table I and II for the 
baryon-baryon-meson vertex and meson interaction vertex, respectively. 



Table I 



Table II 
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The zeroth-order Green's functions of nucleons, deltas as well as mesons corresponding 
to the free Lagrangian density of Eq. (14) are summarized in Appendix A, where the 
distribution functions of negative-energy states are neglected for fermions. They are kept 
for bosons. Considering that we will derive a transport equation for the pion in which we 
only treat the real pion with positive-energy states, we rewrite the zeroth-order Green's 
functions of the pion as 

^ *) = w^[±T, ~ $T) * " ^ ''<*• *>• < 20 » 

A°+-(x, k) = --— 5[k - u)(k)] [1 + Ux, k)} , (21) 
A°-+(x, k) = -— 5[k - co(k)} U(x, k), (22) 
here u(k) is the energy of the pion. 

III. DERIVATION OF THE QUANTUM TRANSPORT EQUATION FOR 

PIONS 

A. Dyson Equation and Pion Self-energy 

With the discussions of Sect. II we can write down the pion Green's function in the 
interaction picture as 

iA 7I (l,2)S ij = (T[exp(-i I dxff J (a;))7rj(l)irj(2)]) 

- (T[exp(-i I dxHj{x))iz !{!)}) (T[exp{-i f dxHj(x)) 77^2)}}, (23) 

here i, 2, 3 represent the isospin indices of pion. In the following we suppress this 
subscript because at the end we will obtain a RVUU-type transport equation which is 
averaged on the isospin. Furthermore, the second term on the right-hand side of Eq. (23) 
vanishes in the spin- and isospin-saturated system. By expanding Eq. (23) perturbatively 
one can obtain the Dyson equation for the pion Green's function, which reads as 

iA^l, 2) = (1, 2)+ f dx 3 f dx 4 A£(l, 4)11(4, 3)*A W (3, 2), (24) 

here 11(4, 3) is the self-energy of the pion, which is also a matrix on the time contour 
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, n— (4,3) n-+(4,3) , 

n 4,3) = K J • (25) 

n+-(4,3) n++(4,3) 



Eq. (24) is coupled to the Dyson equation of the nucleon P8| , |52" | 

iG(l, 2) = zG°(l, 2)+ / dx 3 f dx 4 G°(l, 4)£(4, 3)iG(3, 2) (26) 

and delta |5(| 

iG aP (l,2) = zG%(l,2)+l dx,f rfx 4 GL(l,4)S^(4,3)^(3,2) (27) 

through the self-energy terms of 11(4,3), £(4, 3) and £^(4,3). Here G(l,2), G a/3 (1,2) 
are Green's functions of the nucleon and delta, and £(4,3), £^(4,3) are the respective 
self-energies. Eqs. (24), (26) and (27) are a set of dynamical equations for the hadronic 
matter. From Eqs. (26) and (27) we have derived the RVUU-type transport equations 
for the nucleon [pT5|, ^ - p9| , |52| and delta distribution functions. In this work we will 



develop a transport equation for the pion distribution function from Eq. (24), in which 
both the mean field and collision term will be expressed analytically. Since the lowest- 
order Feynman diagrams contributing to the two-body scattering cross sections are the 
Born diagrams, we consider the pion self-energy 11(4, 3) up to the Born approximation. 
In principle, one should calculate the in-medium cross sections and meson fields for all 
the particles within a relativistic G-matrix theory. However, since we have to deal with 
many reaction channels and many degrees of freedom, such calculations seem to be out 
of the present practical possibilities. For a qualitative insight in the cross sections and 
potentials we think that the Born approximation will be sufficient. A comparison between 
the cross sections for a* NN ^ NN and o~* NN ^ NA calculated in G-matrix theory p4] , |65| and in 



Born approximation [47|,4q| shows differences only in the order of 10-20%. 



The pion self-energy up to the Born term can be written as 

n(4, 3) = U HF (A, 3) + n Born (4, 3), (28) 

here n#_p(4, 3) is the Hartree-Fock self-energy of the pion and n£ orn (4, 3) is the Born 
self-energy. The corresponding Feynman diagrams are given in Fig. 2, 3 and 4. 
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Fig. 2 




Fig. 3 




Fig. 4 



In Fig. 3 we only take the baryon-loops into account since the contributions of meson- 



loops (ct-tt and p-ir mixed loop) are negligible [66] at zero temperature (finite temperatures 



are not taken into account explicitly in the present framework of microscopic transport 
theory). Furthermore, since the pseudovector form is chosen for the wNN and 7rAA 
coupling, as discussed in Ref. |5(J (Sec. 8.3) the contribution of the sigma-pion coupling 



term to the ttN S-wave scattering lengths is small of order m^/(M^m a ) and can be 
neglected, we drop the contribution of Fig. 2(a) and Fig. 2(b) to the pion self-energy. 
Therefore, only Fig. 2(c) contributes to the Hartree-term of the pion self-energy, which 
plays a role in the case that a large amount of pions are produced in relativistic heavy-ion 
collisions at very high energy. For the Born term we consider the Feynman diagrams 
contributing to the tt + N ^ tt + N, 7r + A ^ 7r + A and tt + tt — > tt + tt elastic scattering 
processes as depicted in Fig. 4. For the same reason we neglect the contribution of the 
cr-exchange in Fig. 4(a) and 4(b). The Hartree-Fock self-energy 11^(4,3) and Born self- 
energy U.Bom(4:, 3) can then be expressed as 

rw(4, 3) = u h (a, 3) + rr, oop (4, 3), (29) 

IW4, 3) = IWi (4, 3) + II AA -i (4, 3) + IWi (4, 3) + U NA -i (4, 3), (30) 

n Born (4, 3) = n a (4, 3) + n 6 (4, 3) + n c (4, 3) + 11,(4, 3), (31) 

where the lower subscripts N" 1 and A -1 in Eq. (30) denote the particles described by the 
nucleon and delta distribution functions rather than the anti-particles which are not taken 
into account in this work because the gap between the effective masses of particle and 
anti-particle is much larger than the pion mass even at three times normal density (see 
Fig. 5). For even higher densities and temperatures the production of particle and anti- 
particle pairs becomes more important and the anti-particle degree of freedom should 
be taken into account. The lower subscripts a, b, c, d in Eq. (31) denote the terms 
contributed from Fig. 4(a)-(d), respectively. I1aa- 1 (4, 3) (corresponding to Fig. 3(c) and 
3(d)) and 11^-1(4,3) (corresponding to Fig. 3(g) and 3(h)) are usually neglected in 
the investigation of the influence of the in-medium pion dispersion relation on the pion 
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dynamics in relativistic heavy- ion collisions |16|]I^j2l]f22]j26|| . However, it was recently 
reported that more than 30% of nucleons are excited to the resonance-states, especially 
A-resonance, in the Au + Au collisions at an incident energy of 2 GeV/nucleon [|67[| . 
That means that the contributions of Fig. 3(c), 3(d) and 3(g), 3(h) should be taken 
into account once the problem of the in-medium pion dispersion relation is concerned 
in relativistic heavy-ion collisions. To our knowledge, up to now no investigation of this 
have been made in transport theories. Here we note that this effect has been addressed in 
some non-relativistic calculations of the pion self-energy in static nuclear matter at finite 



temperature 68 



The concrete expressions of self-energies in Eqs. (29) - (31) can be written down 
according to the standard Feynman rules 



11^(4,3) = -(g CT7r m 



n 



NN 



~>NN 



5(3,4)/dx' 3 A°(3',3')zA°(3',4) 
^75^(3,4) f 75 G°(4,3)' 



Hr 



n A A-(4,3) = -hi{g AA )Hr [/> 75 G^(3,4) P l5 G°^(A,3) 
n AJV - l (4,3) = |z(g^ 



Hr 



G°(3,4)P fM P u G° (4,3) 



n 



ATA" 1 



(4,3) 



4 



Hr 



P»P»Gl (3,4)G°(4,3) 



(32) 
(33) 
(34) 
(35) 
(36) 



H (4,3) = £ / dxj dx e (r | g^V; | r 4 )A°(4, 3)(r 4 | g^T* | r) 

tr{(t % | g% N T N p | i 5 )G°(5,6)(i 5 | g% N T N p | t 6 )G°(6,5)} 
A°(4,6)Aj(5,3) J D pJ D p , 
H 6 (4, 3) = £ fdxj dx 6 (r | g^ | r 4 )A°(4, 3)(r 4 | g^ \ r) 

tr{(T e | g AA Tf I T 5 )G°^(5,6)(T 5 | g AA Tf \ T 6 )G° pa (6, 5)} 
A°(4,6)Aj(5,3) J D pJ D p , 
H c (4,3) = - £ fdxif dx 6 (r \ g^Y\ | r 4 )A°(4, 3)(r 4 | g^Y\ \ r) 

(r 6 | gtT\ I r 5 )A°(5,6)(r 5 | g^ | r 6 )A°(6,5) 
A° A (i,6)A A (5,3)D A D A , 



(37) 



(38) 



(39) 
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II d (4,3) = - J2 f dxj dx 6 (r | gtT\ I r 4 )A°(4,5)(r 4 | gj^I* | r 5 > 
A°(5,6)(r 5 | g tT\ | r 6 )A°(6,3)(r 6 | g^T B \ r) 

A° A (4,6)A B (5,3)D A D B , (40) 

In Eqs. (37) - (40) A, B = a, p; r, r 4 , r 5 , r 6 represent the isospin of pions, t§, te denote 
the isospin of nucleons and T 5 , T 6 of deltas. The definition of the symbols is given in 
Table I and II. The transformed four-momentum P M in Eqs. (33) - (36) stems from the 
derivative coupling of baryon-baryon-pion vertex used in our calculations. 

B. Kadanoff-Baym Equation 

Introducing the differential operator of the Klein-Gordon field 

A 01 = <9X + ml (41) 



and applying it to the both sides of Eq. (24), with the help of relation [69 



A 01 A°(l,2) = -5(1,2), (42) 

we obtain 

Aox 1 z-A^l, 2) = -i<5(l, 2)- f dx 3 n(l, 3)^(3, 2). (43) 

It has been shown in Sect. II that only two components of A 7r (l, 2) are independent, from 
which the dynamical equations for the distribution function and the spectral function can 
be constituted |70| . Since we will use the quasi-particle approximation in the derivation, 
the spectral function turns out to be a S function on the mass shell. Thus, in the present 
work it will be sufficient to consider only one component of A 7r (l, 2), i.e., A~ + (l, 2), which 
is directly related to the single-particle density matrix in the case of ti = t 2 [§7]]. The 
equation of motion for A~ + (l,2) can be extracted from Eq. (43). Before doing it, let us 
firstly look at the Feynman diagrams in Figs. 2, 3 and 4 which will be taken into account 
under Born approximation. 
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As is well known, the RVUU-type transport equation contains two important ingre- 
dients, i.e., the transport part related to the real part of the pion self-energy and the 
collision term corresponding to the imaginary part. The Hartree-term of Fig. 2(c) only 
contributes to the real part. However, the loop diagrams in Fig. 3 and the Born diagrams 
in Fig. 4 include both real and imaginary part. It should be pointed out that the baryon- 
lines in Fig. 3 denoted by the symbols N or A represent virtual baryons (nucleon or delta) 
when one calculates the real part of the self-energies. They are not on-shell particles. The 
corresponding terms for the a and u self-energies are neglected in the derivation of Ref. 
|55| for the o and uj transport equations because they used the restriction that all Green's 
functions in the Feynman diagrams should be on the mass-shell. Consequently, mesons 
(a and uj) became free particles in their framework. To our opinion, in computing the 
real part of self-energies which mainly relates to the virtual processes, it is not necessary 
to keep all particles on the mass-shell which will essentially give the imaginary part. It 
is well known that the particle-hole excitation is very important for the in-medium pion 
dispersion relation which will certainly have influence on the pion spectra and pion flow 
in relativistic heavy-ion collisions and should be taken into account in any realistic 
transport models for pions. For the imaginary part of the self-energies from one-loop 
diagrams in Fig. 3, we include only the contributions of Fig. 3(e) and (h), which con- 
tribute to the important A-formation process of iV + 7r — > A and A-decay process of 
A — > N + 7r, respectively. The reason is as follows: the contributions of the imaginary 
part of Fig. 3(a)-(d), in which both the baryon-lines are on the mass-shell, correspond 
to the process that a nucleon (delta) decays into a nucleon (delta) and a pion, which is 
forbidden due to energy-momentum conservation (here we don't consider the Cherenkov 



radiation discussed in Ref. [71], this process might be possible at high densities where the 
pion has a space-like four momentum due to its large potential); Fig. 3(f) and 3(g) don't 
correspond to realistic physical processes when the pion has a positive energy. Since the 
matrix elements are the same for both A-formation and A-decay process, we only need 
to calculate the imaginary part of Fig. 3(e) explicitly. In view of the Born diagrams we 
take only the imaginary parts into account and drop all the real parts which in principle 
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are the corrections to the real part of the Hartree-Fock self-energies. The imaginary part 
of the self-energies can be expressed by 11^(1, 3), which is defined as 



n^(i,3) = n^ rn (i,3) + ng(i,3) 



(44) 



here 11^(1,3) represents the imaginary part of Fig. 3(e). The equation of motion for 
A~ + (l,2) can then be written as 



*Ar + (i,2) 



-£dx 3 (i2enr o ; p (l,3))iA;+(3,2) 



to 



dz 3 n^(l,3)t(iteAj + (3,2)) 



to 

t2 



n+^(l,3)-n-+(l,3)l *A;+(3,2 



+ / rfx 3 n-+(l,3) zA+-(3,2)-^A;+(3,2; 

Jt 1 



(45) 



Eq. (45) is the so-called Kadanoff-Baym equation f72fl . Here the symbol Re denotes 
the real part of the corresponding self-energies. The second term on the right-hand side 
of Eq. (45) corresponds to the spreading width in the spectral function. It should be 



dropped under the quasi-particle approximation [|7(| which will be introduced later. The 
structure of the third and fourth term on the right-hand side of Eq. (45) implies that 
they contribute to the collision term of the transport equation. The concrete expressions 
of the self-energies read as 



n 



AA 



.(1,3) 



n^_i(i,3) 



n 



AN' 



41,3) 



-Ag^ a f{ f ^ 3 A°-(3',3')iA°-(3',l) 
4 Jto 

+ j\x' 3 A ++(3',3')iA° cr + -(3' 1 l)}, 

/^G°-(3,l) P l5 G°-(l,3) 
P l5 G°-~ (3,1) P l5 G°-~^(l,3) 
G°-(3,1)P^G^-(1,3) 



MSNNYtr 



-^(g AA ) 2 tr 



n^A-x(i,3) 



■>AN, 



-tr 



(g A7V ) 2 tr [P»P»Gl--(3,l)G*-(l,3) 



ng)(l,3) = ^(g\ N ftr [G ^(3,l)P»P»Glf (1,3) 
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(46) 
(47) 
(48) 
(49) 
(50) 
(51) 



The expressions of the Born terms are rather complicated. If one does not write out the 
isospin factor explicitly, IT^ T (1,3) can be expressed as 



n^(i,3) 



dx 5 f dx 6 A°^(l,3) 

to J to 

{tr [G° tt (5, 6)G 0tt (6, 5)1 A^ ±T (1, 6)A^ T (5, 3) 



1,6)A°^(5,3) 



+tr [G 0±± (5, 6)G 0±± (6, 5)" 
-tr [G° t± (5,6)G 0±t (6,5)" 
-tr [G°^(5, 6)G°^(6, 5)] A°^(l, 6)A°^(5, 3)}. 



A° ±± (1,6)A^(5,3) 



(52) 



Other Born terms can be written down in the same way. 

C. RVUU Equation of the Pion 

Defining X = \{x\ + X2), y = x\ — £2, x' = £3 — X2 and taking the Wigner transfor- 
mation on the both side of Eq. (45), we arrive at 

~\d%d x - iPttf -P 2 + ml + U H (X) + Reu^-jx, P) - l -d x n H (x)tf 
+ foRen^(X,P)tf - ^M^(I,P)<] ^A- + (X,P) 
dye lPy [ T dx' [ui oU (y - x' , X)iA; + (x', X) 

J — CO L 

-U;+(y-x',X)iA:-(x',X) 



(53) 



Here we have adopted the semi-classical approximation, in which the Green's functions 
and self-energies are assumed to be peaked around the relative coordinate and smoothly 
changing with the center-of-mass coordinate. The details of the Wigner transformation 
are given in Appendix B. The Hermitian conjugate equation of Eq. (53) reads as 

-d*d x + %P% -P 2 + ml + U H (X) + Rell^X, P) + l -d x U H (X)d p 



- -d£ReTlT- p (X,P)d* + -d x ReUi-- p (X,P)d^ 



*a; + (x,p) 



//"CO 
dye** J dx' [n^(y-x',X)<A+-(x',X) 



-Ut M (y-x',X)tA- + (x',X) 



(54) 
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We drop the term of d x d* in Eqs. (53) and (54) since it may be viewed as of higher order 
than the other terms within the gradient expansion used in the Wigner transformation. 
If one would keep this term, the pion Green's function could be non-zero for off-shell 
four- momenta JftJ. In this paper we only consider real on-shell pions. The summation of 
Eqs. (53) and (54) gives 



P 2 -ml- U H (X) - ReUT- p (X, P)\ *A; + (X, P) = 



(55) 



and the subtraction of them yields 
{P»d* + \d x U H {X)dl + \d x ReUT- p {X,P)dr - ^ P ReIl k) - p (X,P)d^A^(X,P) 



n+^(x, P)A~ + (x, P) - n-+(x, P)A+-(x, p) 



(56) 



Now we introduce the quasi-particle approximation and dress the masses and momenta 
in the zeroth-order Green's functions appearing in the self-energies with the effective 
masses and momenta. The canonical variables X, P are then transformed to the kinetic 
variables x, p which will be used in the RVUU code for the simulation of relativistic 
heavy- ion collisions. Since the pion is a pseudoscalar particle, we have P M = p^. Medium 
effects are included through the effective mass which is defined as 



™* ( X 'P) = m ir + Rh(x) + Pell, (x,p) 



(57) 



The on-shell condition is guaranteed by Eq. (55) 



p 2 -LU* 2 ( P )=0, 



(58) 



here 



u*{p) = p +ml (x,p) 



(59) 



We further define a distribution function 



iA w + (x,p) 



71 



co*(p) 



p°-uj*(p) f w (x,p,r), 



(60) 



where 
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-7-1-1 1 9ReU loop( X ,P) | 

Zb ~ 1 " 2^{p} dio lpo =^ {p) ■ (61) 
In the following we drop the derivative term in Eq. (61) since it will cause significant 
difficulty in deriving the collision term. In the nuclear medium the quantum numbers of 
the pion can be either transported as a physical pion or as a delta-hole bound state. In 
several studies of the non-relativistic delta-hole model one considers the mixing between 



these two branches of the pion dispersion relation |21],[E|[26[] . In this case strength is 
redistributed between the two branches as a function of momentum. Therefore, the wave 
function renormalization factor is essential and in principle it can be calculated from the 
energy dependence of the pion self-energy. However, in practical application in transport 
theories only phenomenological simulations of this mixing have been investigated, Since 
we neglect the delta-hole branch in our relativistic dynamical treatment, we put Zg 1 = 1. 
This does not cause any difficulty with the conservation laws as can be seen from Eq. 
(F3). We simply obtain a different but still conserved current. An improvement over 
this not very satisfactory situation might be achieved if one studies relativistic transport 
theories beyond the quasi-particle approximation ||70|| . Especially, the inclusion of bound 
states in transport theory is, however, studied only in a very few non-relativistic case near 



equilibrium up to now, e.g., for the formation of deutron in nuclear matter |74|. It is clear 
that the relativistic bound-state problem is much more involved than the non-relativistic 
one. Therefore, we neglect this problem here. Through inserting Eq. (60) into Eq. (56) 
one obtains the self-consistent RVUU equation for the pion distribution function 

= F c (x,p). (62) 

In Appendix F we show that this equation satisfies the conservation laws of current and 
energy-momentum tensor. The left-hand side of Eq. (62) is the transport part and the 
right-hand side is the collision term, which includes two parts, 

F c (x,p) = F Nw ^ A (x,p) + F d (x,p) (63) 
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stemming from the Nix — > A process and 7r-hadron elastic scattering processes, respec- 
tively. Other reactions are not included in the present work. The collision term can be 
further expressed via in- medium differential cross sections (Sect. V). 



IV. CALCULATION OF THE MEAN FIELD 

In Sect. Ill we derived the RVUU-type transport equation for the pion distribution 
function. The left-hand side of the equation is the transport part and the right-hand 
side is the collision term. The heart of the equation is the mean field, which relates to 
the in-medium pion dispersion relation, and the 7r-relevant in-medium differential cross 
sections. In this section and the next section we will evaluate the concrete expressions 
of them. Before coming to it, we would like to emphasize again that in the present work 
we consider only the n meson as a real meson, a, cu and p mesons are still viewed as 
virtual ones. In other words, the terms relating to the distribution functions of a, uo and 
p mesons vanish. After Wigner transformation, Eqs. (46) - (50) turn out to be 



U H (x 

n 
n 



° )2 I ^ A ^~(^)* A °"M), 



AAj 



NN- 1 V X 5 P) 

x,p) 

4 

4 



(2^) 4 
2 /* A 



d 4 q 



tr 



(27T) 4 

d A q 



tr 



tr 



n 



NA- 1 



(2»)- 



fry 5 G° (x, q) fry 5 G° (x, p + q) 
^ 5 G°--(x,q) ^G'-'^p + q) 
G°-(x,q)p^G ~~(x,p + q) 
P tl p U Gt-(x,q)G -(x,p + q) 



(64) 
(65) 
(66) 
(67) 
(68) 



In the next step we insert the zeroth-order Green's functions for baryons (Appendix A) 
and pion (Eqs. (20)- (22)) into Eqs. (64) - (68) to obtain concrete expressions of the real 
part of the pion self-energies. Several approximations are made here. Firstly, we take the 
quasi-particle approximation in which the free masses and momenta in the zeroth-order 
Green's functions are addressed by the effective masses and momenta. Secondly, the first 
term on the right-hand side of Eq. (20), which will appear in Eq. (64) in the calculation 
of the Hartree term, is dropped as usually done according to the physical arguments 
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(otherwise, it will cause divergence) |75|. Thirdly, in computing Eqs. (65) - (68), we drop 
the contributions of anti-particles contained in the baryon Green's functions of G° (x, q) 
and G°~~(x,g). The zeroth-order Green's functions used in this section then read as 

1 



G°-(x,q) 



q l — m l a + ie ' 



u*(q) 
q 1 +m* 
2E*(q) 
<j+m A 
2E* A (q) 



10 



w*(g)]/,r(x,q,r 
1 



g - E*(q) + ie 



+ 2iri8[q -E*(q)]f(x,ci,T) 



DM 



1 



go - E* A (q)+ie 



2ni S [g - E A (q)] / A (x, q, r) 



(69) 
(70) 
(71) 
(72) 



where E*(q) = [q 2 + m* 2 } 1 ^ 2 , E A (q) = [q : 



.21 1/2 



Dfj, v (q) is given in Appendix A. The 



definition of m* and m* A will be given later. It is interesting to notice that only the 
Green's functions on the upper branch of the time contour, which are similar to the ones 
used in the standard effective field theory, enter in the calculations. 
The Hartree term can be directly worked out 

1 



n H {x) = 

here ps( 7r ) is the scalar density of pion 

( \ 3 



r) Ps(7r), 



(73) 



dq 



=/ w (x,q,r). 



(74) 



For the one-loop diagrams we have to distinguish the real and virtual baryons. The first 
terms on the right-hand side of Eqs. (71) and (72) describe the virtual nucleon and delta, 
which are denoted by the N and A on the Feynman diagrams in Fig. 3. The second terms 
with distribution functions represent the real nucleon and delta, and denoted by the iV -1 
and A -1 on the Feynman diagrams. Through inserting Eqs. (71) and (72) into Eqs. (65) 
- (68), after some straightforward algebra we obtain the real part of the self-energies 



i?en--_ 1 (x,p) = -2(g^) 2 



Rell-.^p) = -5(gl A ) 2 J ■ 



d 3 q 
(2tt) 3 
d 3 q 



A(p, q) 



/(x,qr) 



2tt 



AE*(q)E*(p + q) E*(p + q) - E*(q) - p 
B(p,q) / A (x,q,r) 



+ Po 



AE* A (q)E A (p + q) E* A (p + q) - E A (q) - p 



+ Po 



-po 



-Po 
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ReU AN ^(x,p) 
ReU~ A ^(x,p) 
here 



r>ANj 



d 3 q 

d 3 q 
(2^)3 



C(p, q) 



/(x, q, r) 



4E*(q)E* A (p + q) E* A (p + q) - E*(q) - p 

D(p,q) /a(x, q, r) 

4E A (q)E*(p + q) E*{p + q) - E* A {q) - p 



Po 



+ Po -> - 



Po 
Po 



A(p, q) 
(P-Q)a 
B(p, q) 

(P-Q)b 
C(p,q) 

(p-q)c 

D(p,q) 
(P-Q)d 



pI(P-Q)a - 2m* 2 pl + 2(P-q) A 



(E*(q)+p ) 2 -E* 2 (p + q)-p 



i 2 Pu m *A (P-q)B-m* A -2p Jp.q) B (p.q) B -3m A 



-bplmt [(p.q) B ~ 2mt\ - 2(p.g)| [2(p.g)| + 2(p.q) B m* A 2 + 5m*^\}, 
1 
2 



(^(g)+Pbr-J5S'(p + 9)-p{ 



3mf 



+ 2 pI(P-Q)c - m A p 2 + (p.g)c (p.q) c + m* 2 + m*m^ 



(E*(q)+p ) 2 -E A 2 (p + q)+m* 2 



m 



*2 



Pi 



3m*? 



m *APu~ (P-Q)d (p-q)D + m*m* A +m A ' 



(E* A (q) + p ) 2 - E* 2 (p + q) + m* 2 - m*£ - p 



(79) 
(80) 

(81) 
(82) 

(83) 

(84) 
(85) 

(86) 



Here we already dropped the contributions from virtual particle-particle excitations 
(which are divergent) in consistently with the mean field approximation. Otherwise, 
one has to renormalize it which may be difficult in many situations. The effective mass 
of pion is defined in Eq. (57). The effective masses and momenta of nucleon and delta 



are defined as [4S . 3Q, 



m*(x) = M N - g a NN a(x), 
m* A {x) = M A - g AA a(x). 

r a — 



If) 



Pa( x 



(89) 

(90) 



The mean fields of a(x) and u^{x) are obtained through the following field equations 
within the local density approximation 
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m 2 a a{x) + b(g a NN ) 3 a 2 (x) + c(g^)V 3 (x) = g a NN p s (N) + g£ A Ps(A) + -g a7r m a p s (7r), (91) 

ml^(x) + {g% f 2 ml {^{x)f = g% N p» v (N) + g^(A). (92) 

The scalar and vector densities of the nucleon, and delta are defined as 

PS® = M~J ^ , ^ , /<(*, q, t), (93) 



(2tt) 



+ m* 



7W /" , <? M 



PvH) = 777-^ / c?q =/i(x, q, r). (94) 

( 27r ) J ^Jq 2 + mf 

The abbreviations i=N, A, and 7(i)= 4, 16, correspond to the nucleon and delta, respec- 
tively. 

In Eqs. (87) - (90) we have dropped the Fock term of the nucleon and delta self- 
energies since it makes only a small contribution. The Feynman diagrams for the Hartree 
term of nucleon and delta self-energies can be drawn in the same way as in Fig. 2 by 
replacing the pion external line with the nucleon and delta line. One may notice that 
the main contributions to the mean field of the nucleon and delta come from the Hartree 
term while to the pion from the Fock term (one-loop diagrams). The situation is caused 
by the pseudovector coupling for the pion adopted in our considerations. If one uses 
pseudoscalar coupling for ttNN and 7rAA vertex, the pion will have a scalar self-energy 
from the Hartree term similar to the nucleon's and delta's. But this term turns out to be 
so large that the effective mass of the pion will become almost five times more massive 



at the normal density of nuclear matter than in the vacuum ||66fl ! This may not be the 
realistic case. 

V. CALCULATION OF THE COLLISION TERM 

A. In-medium Nn — > A Cross Section and A-decay Width 

Now we come to calculate the right-hand side of Eq. (62), i.e., the collision term. The 
part corresponding to the A-formation cross section reads as 
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Ntt — >A 



[x,p) 



U t(e)( X iP) A n + ( X 'P) -n 3( 5(a;,p)A+ (X,p) 

d 3 q f d 3 k 



(95) 



(2nfJ (2tt) 3 

here F®, F® are the Nordheim-Uehling-Uhlenbeck factors of the gain (F 2 °) and loss (Fj 3 



terms, 



F 2 ° = [1 + ^(x, p, r)] [1 - /(x, q, r)] / A (x, k, r) 
-^1° = A(x, P, r)/(x, q, r) [1 - / A (x, k, r)] . 



g, fc) is the transition probability of the Nn —>■ A process 



(96) 
(97) 



W(p,q,k) 



[San) 



QLU*(p)E*(q)E* A (k) 
(San? 



tr[^+m*)p^p v (^ + ml)D^(k)) 



1 



(m A + m ) -m n ( 



* \2 *2 

,m — m A ) — 



(9£ 



uJ*{p)E*(q)E* A (k) l%mf 

In working out the second equality of Eq. (98) we have used the relations k = p + q 
and p.q = (m A — m* 2 — m* 2 )/2 from the energy-momentum conservation and on-shell 
conditions. 

In the above derivation all baryons are treated as elementary particles as usually 
done in quantum field theory. However the delta is a resonance that can decay. A mass 
distribution function of the Breit-Wigner form is commonly introduced to describe the 



resonances with broad-widths p9|j76[] . However, one mostly discusses the problem in the 
free space. Here we assume that the same form of distribution function applies to the 
medium with the free quantities replaced by the effective quantities. That is, we introduce 
a Breit-Wigner function for the A-resonance in the medium 



F(mt) 



m*T*(\ p |) 



(99) 



7r(mf-mf) 2 + mfT* 2 (\ p |)' 

where T*(| p |) is the in-medium momentum-dependent A-decay width, p is the relative 
momentum between nucleon and pion in the A-rest system 



[m A — (m* + ml) 2 ] [m A — (m* — m*) 2 ] 
4m A 2 



(100) 
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m,Q is defined by Eq. (88) with the free A-mass M A replaced by its resonance-mass Mo. 
Inserting the mass-distribution function of Eq. (99) into Eq. (95), we have 

d 3 q r d 3 k 



F N ^ A (x,p) = J^l^j dmZ(2n)^\p + q -k)W(p, q ,k)F(m%)(F°-F°) 

(101) 



(27r) 3 i (2tt) 3 
' rq ;va abs (s)(F°-F°). 



In the second line of the above equation we already expressed the collision term with the 
cross section [77]]. Since we are now in the A-rest system, the effective total energy of the 
system s equals to the effective A-mass m* A . cr afes (m A ) reads as 



crab Am, 



2tt (g? 2 



iAN) 



(m* A + m*y - ml 2 (m* — m* A ) — m* F(m* A ). (102) 



| p | 9m A 

Performing an average over the initial states and writing out F(m* A ) explicitly, we arrive 
at the cross section of the Nit — > A process 



a 



NTT^A{ m A 



Ian) 2 



18m*^ | p 



f * i *\2 *2 

[m A + m ) -TO,; 



2 r 



( * * \2 *2 

\m — m A ) — 



x- 



m* T*(\ P I) 



mf - mf) 2 + mfY* 2 {\ p |)' 



(103) 



According to Ref. IjTqj , the resonant cross section can also be expressed by means of the 
decay width 

4tt (2J+1) (2/+1) m* 2 r* 2 (| p |) 



(104) 



N *~* A | p | 2 (2s + 1) (2t + l)(2r + 1) (mf - m* 2 ) 2 + m£ 2 r* 2 (| p |) ' 

here J, J are the isospin and spin of the A, t, s are that of the nucleon, and r is the 
isospin of the pion. Comparing Eqs. (103) and (104) we obtain the in-medium A-decay 
width, which reads as 



r*(| P 



IPl (San) 2 
16nmQ 6m A 



( * i *\2 *2 " / * * \2 *2 

(m A + m ) — m n [m — m A ) — m n 



(105) 



In this way, we give a clear-cut relation between in-medium A-decay width and A- 
formation cross section. 



B. Elastic Pion-hadron Scattering 
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In this subsection we derive the analytical expressions for calculating the in-medium 
tt + N ^ n + N, 7r + A^7r + A and it + n — > n + ir elastic scattering cross sections. The 
corresponding part of the collision term can be written as 

Fei{x, p) = - [u+- rn {x,p)A- + {x,p) - Ust n (x,p)A+-(x,p) 
r d 3 p 2 r d 3 p 3 r d 3 p A 4 (4) 

xW(p,p 2 ,p 3 ,p 4 )(F 2 -F 1 ), (106) 

here F 2) F 1 are again the Nordheim-Uehling-Uhlenbeck factors 

F 2 = [l + / w (x,p,r)] [l±/ff 2 (x,p2,r)]/ff 3 (x,p3,r)/ff 4 (x,p4,r), (107) 
F 1 = / w (x, p, r)/fl- 2 (x, p 2 , r) [1 ± /h 3 (x, p 3 , r)] [1 ± /h 4 (x, p 4 , r)] , (108) 

if 2 , #3, #4 can be 7r, N or A; the symbol + assigns to bosons and — to fermions. The 
transition probability W(p,p 2 ,p 3 ,p4) for different channels reads as 

w^ M , n , Pi) = lfo . (p) ^i; )£ , (P4) ^. (no, 

W^^{p,p 2 ,pz,pi) = V M ) 4 T C $ C + (g^g* ) 2 T d $ d 

lQuj*(p)uJ*(p 2 )uj*(p 3 )uj*(p 4 ) ^ L J 

+P3^P4- (HI) 

where A, 5 = a, p. T a ^ d is the isospin matrix and $ a -d is the spin matrix. The 
subscripts a,b,c,d denote the terms contributed from Fig. 4(a)-(d), respectively. The 
concrete expressions for T a _ d and <& a -d are 
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= E ( r 1 
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(112) 


T b 


= E <r 

T 2 r 3 T 4 


1 

1 'p 
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k)(T 4 


\t p \t 2 )(t 2 \t p \t 4 )d;d p , 
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T c 


= E 

r2r3T4 
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3> a = (HfD^tr [ 7p (j* 2 + m*) 7p (i4 + m*)] 



$ 6 = (^DWtr [ lp {fo + m A )D°<>(p 2 ) lp (fa + m A )D pa (p 4 )} 



x 



1 



n 2 



(p -P3) 2 - 



(lAJ ^A^A 



1 



(P -P3) 2 - m A 
1 



(116) 

(117) 

(118) 
(119) 



ijP ~ P3) 2 -rn 2 B (p- p 4 ) 2 - m\ ' 
We further express the right-hand side of Eq. (106) by the differential cross sections [77 

d 3 p 2 



Fel(x,p) 



(2k)- 



: v 7T a n (s,t)(F 2 - Fi)d£l, 



(120) 



here a n (s,t) represents the in-medium differential cross sections oi n + N — > 7r + AT, 
7r + A — > 7r + A and 71 + 71 — > 7r + 7r elastic scattering. Its concrete expressions can 
be obtained through computing Eqs. (112) - (119) and finally transforming it into the 
center of mass system. We give the explicit expressions of ct w n->tvn(s, t), o^a^ttA^, t) and 
a WK - >mc (s,t) in Appendix C. After averaging over initial states, the in-medium total cross 
sections can be calculated through the following equations 



ttN — >7r7V 



-J C r wN^TTN(s,t)dQ, 
nA-nrA = Yq / VirA-nrA{s,t)d£l, 



a 



7T7Y — *7T7T 



q I — 7T7T 



[a, t)dQ. 



(121) 
(122) 
(123) 



Of course, in calculating the cr^A->7rA one should also take the A-decay width into account. 
However, the strict treatment of Breit-Wigner distribution function as in Sect. IV might 
cause complexity in the derivation procedure since we are now concerning two deltas in 
a scattering process. Practically, we usually introduce a centroid A-mass in numerical 
calculations which can include the influence of A-decay effectively. For the detailed de- 

J. At the end we can rewrite Eq. (62), 



scription of the method we refer to Refs. |^q,pO|,p 
the RVUU-type transport equation of the pion, in the following form 
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{P^ + \d^ H {x)d^ + ^ReU^p)^ - iff ifen^fo pK} ^^ T) 

= J -^va ah s{s){F° - F°) + J ^v^{s,t){F 2 - F 1 )dSl. (124) 

The first term on the right-hand side of Eq. (124) stems from the Nn — > A process 
(it is angle independent in the center of mass system) and the second term represents 
the 7r-hadron elastic scattering. This is of course interesting to investigate the 7r-hadron 
inelastic scattering processes, which may warrant further studies. 

VI. NUMERICAL RESULTS AND DISCUSSIONS 

In this section we present our numerical results for the in-medium pion dispersion 
relation, A-formation cross section and momentum-dependent A-decay width. The cal- 
culations are performed in symmetric nuclear matter at zero temperature. The baryon 
distribution functions in Eqs. (75) - (78) and (93), (94) are replaced by the corresponding 
9 functions. The coupling strengths of g^v, g% N and b, c, Z are determined by fit- 
ting the known ground-state properties for infinite nuclear matter. In this work we take 
the parameter set 2 of Ref. @, which gives g^ JV =11.77, g^ iV =13.88, 6(g^v) 3 =13.447, 
c(g^ riV ) 4 =10.395 and Z=3.655. The corresponding saturation properties are: binding en- 
ergy E bin =— 15.75 MeV, saturated effective nucleon mass m* N /M N =0.6, compressibility 
if =200 MeV and the saturation density is po=0.1484 fm~ 3 . 

For the coupling strengths of g^ A and g%&, no direct information from experiments 
is available. For simplicity, we employ the argument of universal coupling strengths, i.e., 



Saa = Saw an d gAA = Snn H3I ■ Other choices of the delta coupling strengths now in the 
literatures [0,^1 are no ^ considered here since they give unreasonable results for the pion 
dispersion relation. In this case, we don't have a nonzero delta distribution in relativistic 
mean field calculations [T75| , |5I}| . The contributions of Eqs. (76) and (78) to the real part 



of pion self-energy vanish in the present calculations of the pion dispersion relation. It 
is, of course, not the realistic situation of the dynamical process of energetic heavy-ion 
collisions, where a rather large amount of nucleons are excited to the resonance states |)7j . 
The contributions of Eqs. (76) and (78) will certainly enter the pion dispersion relation 

28 



and might play an important role because of the large spin- and isospin-factor of the 
A-resonance. Therefore, for the use of the in-medium pion dispersion relation presented 
in this work in the study of high-energy heavy-ion collisions, it should be viewed as a 
preliminary step approaching to the realistic description. 

For the coupling strength of glf N , we take the most commonly used value / 2 /47r=0.08 



56fl . The coupling strength of g^ N can be fixed through using the Eq. (105) in the free 
space. If one takes M N =939 MeV, M A =1232 MeV, m 7r =138 MeV and the empirical value 
of r =115 MeV, it turns out /* 2 /47r=0.362, very close to the commonly used value 0.37 
f40|l (if one uses this value in Eq. (105), it gives r =118 MeV, still in the error bar of 
experimental data). In computing the real part of the pion self-energy we use a cut-off 
factor of A 2 = exp(—2p 2 /b 2 ) with b=7 as usually done [|I~6| , |iT . 



A. In-medium Pion Dispersion Relation 

Fig. 5 displays the gap between the effective masses of particles and anti-particles at 
different densities. One can see that the mass gap is much larger than the pion mass even 
at three times normal density. That means that the anti-particles contribute only to the 
high branch of collective modes, which can be neglected in the present consideration. 

Fig. 6 depicts the in-medium pion dispersion relation calculated with Eq. (59) at 
normal density. Here we already drop the contribution of the Hartree term since we do 
not expect a large pion distribution at the energy and density considered. The areas 
of non-vanishing imaginary part (NIP) of the pion self-energy are also indicated in the 
figure. The imaginary part of the pion self-energy in nuclear matter can be derived 
through inserting Eqs. (71) and (72) into Eqs. (65) - (68) and taking both the baryons on 
the mass-shell. The analytical expressions are given in Appendix D. It should be pointed 
out that in calculating the pion dispersion relation of Eq. (59) the self-consistency is 
realized only in the real part of the pion self-energy. From the figure one can find that 
the A-h branch is above the area of NIP. The sound branch is buried in the region of NIP 
contributed from the nucleon-hole excitation. For the pion branch, it firstly increases with 
the increase of momentum and then disappears in the region of NIP contributed from the 
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A-hole excitation when the momentum is larger than three times pion mass. That means 
that at larger momenta the pion can be bounded in the A-resonance. This scenario is 
commonly used in transport models. In the figure one can also find another dotted curve 
buried in the region of NIP from A-hole excitation, which is related to the case that a 
formed-A decays immediately. 



Fig. 5 




Fig. 6 




Fig. 7 



Since we are mostly interested in the pion branch, in the following we discuss it in 
detail. In Fig. 7(a) we show the pion dispersion relation at different densities. The 
corresponding pion self-energy is displayed in Fig. 7(b). At p — 2p one can see the 
numerical unstable since at that density the pion self-energy is very large compared to the 
pion mass. Contrary to the results of the non-relativistic model \ W\ where the pion branch 



always starts from the point of u* = since the pion self-energy is an explicit function of 
p 2 in that model, our results exhibit that the pion dispersion relation (pion branch) has a 
rather different behavior for different momenta. At lower momenta the pion has a positive 
self-energy, which causes the in-medium pion dispersion relation to be harder than the 
free one. The pion self-energy decreases with the increase of momentum. When p exceeds 
the point around 100 MeV the self-energy becomes negative and correspondingly the 
dispersion relation becomes softer than the free one. One may argue that nucleons are not 
very relativistic in nuclear matter, little difference is expected, for slow pions, between the 
non-relativistic and relativistic results. In Appendix E we reduce our relativistic formulae 
to the non-relativistic limit. It is shown that the relativistic effect stemming from the 
Fermi motion of nucleons is negligible. But there does exist an evident difference to the 
non-relativistic model mainly coming from the relativistic kinetics where p 2 instead of p 2 
is used. The dispersion relation of Fig. 7(a) may provide a possible explanation to the 
pion spectrum over whole energy range. In Ref. |26] it was shown that the pion yield is 
overestimated at low momentum whereas underestimated at high momentum when a free 
dispersion relation was used. Fig. 7(c) displays the contributions of different excitation 
modes to the real part of the pion self-energy. One can find the main contribution comes 
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from the A-hole excitation as expected. The self-energy from the nucleon-hole excitation 
always has a positive value whereas the one from the A-hole excitation changes its sign 
from positive to negative at certain momentum point, which controls the behavior of the 
pion dispersion relation. 



Fig. 8 



The above calculations are performed through considering the lowest-order Feynman 
diagrams. The short-range correlations have not been included. In non-relativistic model, 
the short-range correlations are taken into account by means of the Migdal parameter g'. 



This method is frequently employed in the relativistic model [£3] although it has never 
been checked carefully. We follow this way and take the pion self-energy under Random 
Phase Approximation as 

iMI z ~~ {x, p) 

{p " ] = r+(iv*jfeiWM" (125) 

With g'=0.6 we have recalculated the pion dispersion relation which is plotted in Fig. 8. 
The behavior of the pion dispersion relation becomes quite strange and difficult to be 
understood. Similar results were obtained in Ref. [p3j] . That might mean that, it is 
unsuitable to incorporate the short-range correlations in a relativistic model by means of 
a non-relativistic approach. A fully self-consistent inclusion of correlation effects might 
be necessary, which is, however, quite complicated and needs to be discussed in a separate 
paper. Another possibility is that the effective Lagrangian of Eq. (19) might be valid 
only at lower order. 

For a convenient use in the study of heavy-ion collisions, we parameterize the results 
of Fig. 7(a) as 

1.10398 + 0.0790471p + 0.232015p 2 - 0.049101p 3 p = 0.5p 
1.32175 - 0.13706p + 0.165035p 2 - 0.0173563j9 3 p = p ■ (126) 



10 



1.56304 - 0.585238J9 + 0.261195p 2 - 0.0370091p 3 p = 2p 
The unit of uj* and p is m n . 
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B. A-formation Cross Section and A-decay Width 



Now let us turn to the collision term of the pion transport equation. In this subsection 
we study the ttN — > A cross section and the momentum-dependent A-decay width, which 
are the most important channels for pion absorption and production, both in the free space 
and nuclear medium. As has been pointed out in Sect. V, the delta is a physically decay 
particle. A Breit-Wigner function is commonly introduced to describe the broad-width 
of the A-resonance when one considers the A-relevant scattering processes [|1^,[1^,|7(J . 
Consequently, the mass of the A has a distribution with respect to the total energy of the 
system. But in the framework of an effective field theory one only treats a point particle 
with fixed mass. If one introduces an energy- dependent mass, correspondingly, one should 
introduce certain corrections on the interaction vertex. Brueckner [EI suggested a vertex 
function of 

F -{ 1 + ) (127) 

to fit the phase shift of a broad resonance away from the resonance position. Here p n is 
the relative momentum between nucleon and pion in the A-rest system. R is the radius 



of the boundary of the internal region. The typical value for the A is 0.98 fm fT8fl . In 
the same spirit a mixed version of the form factor for the iVA7r vertex was used in our 
previous works [48,|50],52] which has essentially similar effects. More conveniently, here we 
phenomenologically introduce a mass-dependent coupling strength (MDC) g^ N (M/±) = 
g^F, which will be used in the following calculations. The original coupling strength 
g& N is afterwards referred to mass- independent coupling strength (MIC). 

Fig. 9 displays the n + p — > A ++ and ir~p — > A cross section in free space. The dots 
are the experimental data from Ref. [Q. The solid and dashed curves are our numerical 
results calculated with the mass- dependent (solid) and mass-independent (dashed) cou- 
pling strength, respectively. The results with MDC can reproduce the experimental data 
nearly perfectly. Furthermore, our calculations are almost parameter- free. Only g^ was 
fixed by fitting Tq = 115 MeV. That implies that our theoretical framework for describing 
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the pion is quite reasonable although it should be further checked in relativistic heavy-ion 
collisions. 



Fig. 9 Fig. 10 Fig. 11 

In Fig. 10 we depict the momentum- dependent A-decay width in free space, calculated 
with MDC and MIC, respectively. Some phenomenological parameterizations commonly 
used in the transport models are also presented in the figure. Our results with MDC 
are comparable with these parameterizations. But the decay width calculated with MIC 
increases very rapidly with the increase of the pion-momentum, which will open the 
possibility that a A may have a mass much larger than its resonance mass jSPf . This 
may not be the real case. In the following calculations we will use the mass-dependent 
coupling strength for the iVA7r vertex. 

In Fig. 11 we show the in- medium A-formation cross section and A-decay width. The 
effective masses of nucleon and A are determined by Eqs. (87) and (88). The free pion 
mass is used in (a) and (b) while the effective pion mass from the dispersion relation 
of Fig. 7(a) is used in (c) and (d). From the figure one can find the strong medium 
corrections. In the case of free pion mass, the effective A-decay width decreases rapidly 
with the increase of density. The A-formation cross section is enhanced near the A- 
resonance mass but suppressed at other region. The whole shape of curves becomes 
narrower compared to the free one. The in- medium A-decay width was ever studied by 
Kim et al. ]57| where an effective pion mass stemming from the nucleon-hole excitation 
was used. This kind of effective mass is very close to the free mass as stated in their paper, 
and might be seen from Fig. 7(c). They obtained a suppressed A-decay width at normal 
density which is in qualitative agreement with our results. However, in our calculations the 
Pauli blocking of the final nucleon is not taken into account, so a quantitative comparison 
is not possible. The Pauli blocking is of course incorporated in the transport equation 
(124). The effects of Pauli-blo eking alone on the in-medium resonance decay width were 
ever investigated by Effenberger et al. in Ref. |8£ 



In Fig. 11(c) and (d) the medium effects on the pion are incorporated. Compared 
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to Fig. 11(a) and (b) one can find that the effective mass of the pion can change the 
results completely, even the trend of the density dependence. It seems that if the smaller 
effective mass of the pion is taken into account the in-medium A-decay width increases as 
compared to the free mass. The medium effects on the A-formation cross section are now 
exhibited to be important only in the region where the formed-A has the mass around the 
A-resonance mass. It decreases with the increase of density. When the formed-A is far 
away from its resonance-mass, the medium corrections on the A-formation cross section 
are negligible. 

VII. SUMMARY AND OUTLOOK 

A large amount of data with high accuracy for the pion spectrum and pion flow has 
been accumulated. A 7r-beam facility will be available at GSI which will provide new and 
specific data to help to understand the pion dynamics in relativistic heavy-ion collisions. 
However, a self-consistent treatment of pions, together with the nucleons and deltas, is 
still not realized in the current transport models. In view of this fact, we have developed 
a RVUU-type transport equation for the pion distribution function based on an effective 



Lagrangian of the QHD-II model [56]. The closed time-path Green's function technique is 
employed and the semi-classical, quasi-particle and Born approximation are adopted in our 
derivation. We have presented an unified approach to the following problems: First, both 
the mean field (the real part of the pion self-energy) and the collision term (the imaginary 
part of the pion self-energy) of the transport equation are derived simultaneously from 
the same effective Lagrangian and presented analytically. Second, we treat the real pion 
and virtual pion on the same footing. Third, the transport equation of pion is derived 



within the same framework which we applied to the nucleon [r5| , ^7 - |5| , 52 | and delta [5D,pT 
before. Therefore, we obtain a set of coupled equations for the N, A and 7r system which 
describes the hadronic matter in an unified way. 

Within this approach we have investigated the in-medium pion dispersion relation. In 
contrast to the results of the non-relativistic model where a softer dispersion relation over 
the whole momentum range is exhibited, the predicted in-medium dispersion relation 
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turns out to be harder at lower momenta and softer at higher momenta, compared to 
the free one. The main reason for the difference relies on that in our relativistic model 
the pion self-energy has a relativistic kinetics p^ while in the non-relativistic model it 



explicitly depends on three- momentum p ||40|| , in which the real-part of the pion self- 
energy goes to zero when p — > 0. However, a pion in the nuclear medium should in 
principle suffer the interaction of surrounding particles whatever it moves or not. In our 
calculations the real part of the pion self-energy has a positive value at smaller momenta 
f85fl . Consequently, the momentum dependence of the in-medium pion dispersion relation 
becomes very flat and quite different from the free one and that of the non-relativistic 
model. This will certainly have effects on the pion spectrum, pion flow as well as on the 
dilepton production since one of the important channels 7r + 7r~ — > p — > e + e~ explicitly 
dependents on the slope of duj/dp ||27|| . It would be very interesting to check this kind 
of pion dispersion relation in the dynamical processes of relativistic heavy-ion collisions. 
Work on this aspect is in progress. 

Considering that in the nuclear medium the absorption and production of pions are 
mostly realized through the formation and decay of the A-resonance, we have studied the 
A-formation cross section and A-decay width both in the free space and in the medium. 
Our theoretical prediction for the free A-formation cross section is nearly in perfect agree- 
ment with the experimental data. The computed free A-decay width is comparable to 
several phenomenological parameterizations commonly used in the transport models. It 
is found that the effective pion mass has strong influence on the predicted in-medium A- 
formation cross section and A-decay width. It can even change the trend of the density 
dependence. After taking into account the medium corrections on the nucleon, delta and 
pion mass simultaneously, the A-decay width turns out to be enhanced in the medium 
especially at higher momenta, while the A-formation cross section is suppressed around 
the resonance mass. When the formed-A is far away from its resonance mass, the medium 
effects on the A-formation cross section are negligible. 

In investigating the in-medium pion dispersion relation we have taken into account 
only the lowest-order diagrams stemming from the nucleon-hole and A-hole excitation. 
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The effects of short-range correlations have been addressed through following the method 
of the non-relativistic model with the Migdal parameter g'. But the results turn out to be 
quite strange. In fact, the short range correlations can be considered in a relativistic model 
self-consistently by implementing correlation terms like {ip^^sTil)) 2 and {i/jAp,S + 4>) 2 i n the 
effective Lagrangian of Eq. (19). In the mean field approximation one should calculate 
the expectation value of (07^75T0) and (^>A/xS +, 0) which corresponds to the condensate 
of the pion field. Although the detailed investigation of correlation effects may go beyond 
the scope of the RVUU approach, it is nevertheless very interesting to check the effects 
of short-range correlations on the pion dispersion relation even in static nuclear matter 
within a relativistic model self-consistently. It would be extremely interesting if this 
collective instability could be studied in a dynamical situation. 

As usually done in a microscopic transport theory for non-equilibrium system, the tem- 
perature degree of freedom is not taken into account in the present work. One might con- 
sider the temperature degree of freedom in nuclear matter by simply replacing the single- 
particle distribution functions with the Fermi-Dirac (for fermions) and Bose-Einstein (for 
bosons) distribution functions. In this way one can study the effects of temperature on 
the predicted quantities. This is especially meaningful in the present case since the delta 
distribution will enter explicitly at certain temperatures [79|. Then the AA _1 and 7VA -1 
excitation of Fig. 3(c), (d) and (g), (h) might play an important role on the in-medium 
pion dispersion relation. This would be more close to the realistic situation of energetic 
heavy-ion collisions and requires further studies. 

As frequently emphasized in the paper, up to now we only treat the pion as a real 
particle. Other mesons such as a, uj and p still remain to be virtual ones. It is of 
course interesting to develop transport equations for u and p mesons as well as for other 
experimental observable mesons K, K*, 77, ... within the present framework. Among 
them, the p meson is especially important for dilepton production which can not be 
explained by the transport models currently on the market. Medium corrections on the 
properties of u, p as well as 77 mesons may provide a possible explanation p2-34,90-p2|. 

Almost in all practically used RVUU-type transport calculations, the local-density 
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approximation is employed in order to realize the numerical solution of the equation. This 
drops the retardation effects of the mean field, although the equation itself is constructed 
from the relativistic model. From Eq. (2) one can find that it is possible to describe the 
propagation of real as well as virtual a and uj mesons simultaneously. In this case one 
can obtain the time-dependent mean field and the transport equations for the a and uj 
distribution functions at the same time. It would be interesting to study this problem 
and explore it just from a theoretical point of view. 
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In this appendix we present the zeroth-order Green's functions of nucleon, delta and 
mesons used in this work. 
(1) Nucleon 
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APPENDIX A 



(r[^(l)^(2)]>=iG°(l,2)^, 



(Al) 




(A2) 



G° TT (x, k) 




5[k -E(k)]f(x,k) 



(A3) 



G 0+ -(x, k) 



E{k) 



S[k -E(k)} [l-f(x,k)](H + M N ), 



(A4) 



G°- + (x, k) 



E{k) 



S[k -E(k)]f(x, k)(^ + M N ). 



(A5) 



(2) Delta 
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Wam(1Wa„(2)]> = -iG°(l,2)S Tl>T; 



(A6) 



iGL (1,2) = * 



/7° f-r l-V - **^ 1- * 2 ) 
(27r) 4 ^ ' 



G^M) = (# + M A )L> 



/if 



±1 



77/ 



<J[fco-S(A:)]/ A (x, fc) 



Til 



E(k) 



5[k - E{k)} [1 - f A (x, k)](f6 + M A )D flu , 



E(k) 



6{k -E(k)]f A (x,k)(k' + M A )D : 



11 2 



(A7) 
(A8) 
(A9) 
(A10) 
(All) 



(3) Mesons 



(T[$ A (1)$ B (2)]> = iL> A A A (l,2)<^ B , 



(A12) 



where A, 5 = cr, lu, ir, p, 



D A = D^D 1 



A^Ai 



D A and D A are denned and listed in Table I and 



(A13) 



A° A ^(x,k) 



iA° A (l,2)=ij-0^A A (x,k)e 
±1 



k 2 — m\ ± ie 



-ik{x\—X2) 

2m 6 



A° A ±T (x, k) = -2m 5 



k2 - m A\ fA(x,k), 

[e(±ko) + f A (x,k)}. 



(A14) 

(A15) 
(A16) 



Here the number 1, 2 represent xi, x%. ti, t' 2 denote the isospin of nucleons and Ti, T' 2 
denote those of deltas. f(x, k), f A (x, k), f A (:r, k) are nucleon, delta and meson distribution 
functions, respectively. The abbreviation for isospin on the distribution function has been 
suppressed. 

APPENDIX B 

In this appendix we perform the Wigner transformation of Eq. (45), which can be 
easily realized by means of the following formulae [ITS] 
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h(x)g(x,y) 
<?x'f(x,x')g{x',y) 



(-iP» + -d*)f(X,P), 



h(X)g(X,P)- 
h(X)g(X,P) + 



1 3h(X) dgjX, P) 

2 dX^ dP^ 

1 dh(X) dg(X, P) 

2 dX» dP, ' 



f(X,P)g(X,P) + - 

df(X,P)dg(X,P) 
<9A> dP„ 



df(X,P)dg{X,P) 



OP, 



dX» 



here 



X=-(x + y). 

After Wigner transformation the different terms in Eq. (45) turn out to be 

,1 



n H (i)A- + (i,2 



dZd£-iP"dZ -P 2 )A- + (X,P), 



U H (X)A^(X,P) - -d^ x U H (X)d^A- + (X,P), 



(Bl) 
(B2) 

(B3) 
(B4) 



Mr M Jl, 3)A;+(3, 2) — ReU To -(X, P) A;+(X, P) + % - \dpReIlf~(X, P)d* A~ 



n+7 z (i,3)A; + (3,2 



- d^ x ReUr~ p (X,P)d^A- + (X,P) 
Ut J l (X,P)A- + (X,P). 



(B5) 

(B6) 

(B7) 

(B8) 

-(X,P) 
(B9) 
(BIO) 



In Eq. (BIO) we have dropped the contributions from derivative terms. That means that 



collisions are performed at instantaneous time: Boltzmann Ansatz [54 , |72 . 



APPENDIX C 



In this appendix we present analytical expressions of in-medium differential cross sec- 
tions for 7r-hadron elastic scattering. 

(a) Differential cross section of in-medium tcN — > ttN scattering: 

1 



VnN->uN(s, t) 



(27r) 2 s2(t-m 2 p )' 



m* 2 + ml 2 - 2s) 



x (m* 2 + m* 2 ) + s 2 + st- m* 2 t 



(CI) 
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where 

s = (p + P2) 2 = [u m (p) + E*(p 2 )] 2 - (p + p 2 ) 2 , (C2) 
t = \{2m* 2 + 2m; 2 - s) - ^(m* 2 - mf) 2 + 2 | p 1 1 p 3 | cos 9, (C3) 

I P l = l P3 l = 7T~r\l I s ~ m * 2 ~ m * 2 ) 2 ~ 4m* 2 m* 2 , (C4) 

Ay S 

and 9 is the scattering angle in cm. system. 

(b) Differential cross section of in-medium 7rA — > 7rA scattering: 



1 5( g/OT g£ A ) 2 
(27r) 2 s36m*£{t-m 2 ) 2 



^a^a(s, f) = /o ; 9 "'XT**'™ 18ml 4 (mf + mf) 2 - 2mf{l%s + lit) 



-m*£(36m* 2 s + lQmft - 18s 2 - 2Qst - It 2 ) 

+mf(Sm* 2 st - Amft + 2m* 2 t 2 - 4s 2 t - 6st 2 - t 3 ) 

+mft 2 {mf -2s) + st 2 (s + t)}}, (C5) 



where 



s = (p + P2) 2 = [w*(p) + £ A (p 2 )] 2 - (p + P2 ) 2 , (C6) 

t = \(2mt + 2mf - s) - ^(m* 2 - mf) 2 + 2 | p || p 3 | costf, (C7) 

I p | = | p 3 |= — !-= J(s - mf - mf) 2 - 4m* 2 m A 2 , (C8) 
Zy s 

(c) Differential cross section of in-medium nn — > nn scattering: 

a n ^ n (s, t) = [D(s, t) + E{s, t) + {s,t^ u)} , (C9) 



D ^-m^ + w^ (4< - 28 -^ (C10) 



256(t — m 2 )(u — m z a ) 8{t — mj)(u — m l p ) 



( C 2&o-K m a&pir) 



Ami 2 -2s-t t- s 



+ 



(Cll) 



8(t — m 2 ){u — m 2 ) 8(t — m 2 ){u — m 2 ) 

where the function D represents the contribution of the direct term and E is the exchange 
term and 



40 



s = (p + P2) 2 = [w*(p) + ^*(p 2 )f - (p + P2) 2 , 
t = i(4m; 2 -s)+2|p||p 3 |cos#, 

m = 4m* 2 — s — t, 
1 



P3 



v s - 4m?. 



(C12) 
(C13) 
(C14) 
(C15) 



APPENDIX D 

In this appendix we present analytical expressions of the imaginary part of the pion 
self-energy in nuclear matter, 
(a) For space-like p^ 



ImU NN _ 1 (x,p) = 



*2 2 

^(g« NN nE F - E *), 

7T I p I 



(Dl) 



where 



E F = (m* 2 + k\fl\ 



E*=min (E F ,E max ), 

-Ejm.n.T 771 CI X 



1 1 4m* 2 
m*, E F - I Po I, -- I po I +- I p I (1 — ) 1/2 

2 2 p£ 



(D2) 
(D3) 
(D4) 



ImU 



AA- 



5o 2 

(*,p) = 18nlp ] m *M - H< + iomi 4 )(gi A ) 2 (^ - 



(D5) 



and 



mt + {k F ) 2 



nl/2 



A — mm , E, 



A ) 

max/ 1 



Ei 



max 



m* A ,E£- 



Po 



1 

2 



^|P|(1-^) 1/2 
2 P 2 , 



(D6) 
(D7) 
(D8) 



12vr I p 



f C(p,g)(£ F -£*), 



(D9) 



where 
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E* = min (E F ,E max ), 



E m ax TflCLX 



* 7-1A i i 1 i i 1 i i/2 4?7i ^\i/o 

m p ,~o Po e+ o P ( e a") 

2 2 



and 



e = 



p 2 + m* 2 - m*l 



C(p, q) 



2 / * *\ 2 

P» ~ ( m A - m ) 



'2 ( * i *\ 2 
*V - ( m A + m ) 



3mf 

The definition of E F and £^ is the same as in Eqs. (D2) and (D6). 

V 7T \2 



ImU-.^p) = ^fL-C(p,q)(E^ - E*), 

1Z7T | p | 



where 



E* = min (Ep,E max ), 



E m ax TflCLX 



™a, E F - | po |, 4 bo I e> + ± | p | (e' 2 - 



P 2 * 



and 



(b) For time-like p M 
where 



, v\ + m A — m * 2 



{ -^C(p,q)(E u -E d ) 



*\2 



12tt|p| 



< p 2 < (to a - to*) 
otherwise 



E u = rain 
E d — min (Ep, E, t 



1 , 1 , 

e f, -g I Po I e + g I P 



| m *2 
fg 2 _ ^1/2 



p; 



E, mn ™ max 



TO*, | po 



I 1 I I 1 I 1/2 4m * 2 M/2 

l,-2l^|e--|p|( e 2 -- r )V 2 



(D10) 
(Dll) 



(D12) 
(D13) 



(D14) 



(D15) 
(D16) 



(D17) 



(D18) 



(D19) 
(D20) 
(D21) 



For the time-like p M , the contributions of lm\i~^ N _ 1 {x ,p) and IraIl A ~ A _ 1 , which describe 
the particle-antiparticle decay processes, vanish since we neglect the anti-particles in the 
present framework. In this work the numerical calculations are performed in cold nuclear 
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matter with the assumption of chemical equilibrium, i.e., Ep = Ep when a delta is 
produced. In this case the A-decay process is Pauli-blocked since a delta can only decay 
into a pion and a nucleon with a momentum smaller than the nucleon Fermi-momentum, 
which leads to ImU~j v ~ A _ 1 ( y x,p) = when p >| p |. 

APPENDIX E 



In this appendix we introduce the non-relativistic approximation for Eqs. (75) and 
(77). These two terms stemming from the particle-hole and A-hole excitations are com- 
monly considered in the non-relativistic approach |]40|JT6| , ^T|| . The effective masses and 
energies in the expressions are replaced by the corresponding free ones. In order to make 
a complete non-relativistic reduction we have to start from the full Green's functions 
including the contribution of anti-particles. Eqs. (75) and (77) then read as 



i?en--_ 1 (x,p) = 2(g^) 2 | 



ReU 



AN 



-i{x,p) 



■>AN) 



d 3 q 
(2^)3 

d 3 q 
(2^)3 



A(p, q) 



/(x, q, r) 



2E(q) (p + E(q)) 2 -E 2 (p + q) 

C(p, q) /(x,q,r) 

2E(q) (po + E(q)y-Ei(p + q) 



+ Pn 



(El) 
(E2) 



with the A(p,q) and C(p,q) defined by Eqs. (79) and (83). After some algebra we obtain 



ReU 



NN~ 



ReU 



(x,p) 
(x,p) 



9 



4(g 
(gAAr) 



d 3 q M%pI 
v v ' 1 (2vr) 3 E 2 (q) 
i i d 3 q /(x,q,r) 



/(x,q,r) 



(2vr)3 E(q) 



{p.qf-Mlpl | 2plM N (M A + M N ) 



pl/2E(q)-p-q/E(q)+p 
+ 



(E3) 



Ml 



Mi 



(p.q) 2 - M 2 N p 2 



(M a + M, 



N i 



PI 



2MlE(q) p 2 J2E(q) - p ■ q/E(q) - (M| - M 2 N )/2E(q) + p 



+Pv -» ~Pn\ 



(E4) 



which are the same as Eqs. (8) and (11) in Ref. Taking the non-relativistic limit 

E(q) « M N , Eq. (E3) becomes 



ReU 



NN 



-l{x,p) = {g NN fpl 2 2 ^ P 2 pN 
Po ~ 



(E5) 



with Up = p 2 /2Mn- If one further neglects the relativistic kinematics, i.e., p 2 ^ — > — p 2 , it 
returns to the standard non-relativistic formula stemming from the particle-hole excitation 
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40|]16|1 . Therefore, the relativistic effects stay in two aspects: one is the Fermi motion of 
nucleons in a nucleus which is small; another one is the relativistic kinetics which turns 
out to be substantial in our calculations as can be seen later. In the no n- relativistic limit, 
Eq. (E4) becomes 

ptt— / \ 4 , 2 M N p 2 8 2 {Ma + Mn) 2 

RetlAN-ifap) = -{g AN ) -J^PN + g(gA7v) PftPN 

1 n (M A + M N ) 2 -pl 2u r 2 

~9^ n) Ml pf^l P Pn > (E6) 

with 

Pl Mi M{ 

UR -2Mi;~ 2M N ■ {E7) 

The first and second terms on the r.h.s. of Eq. (E6) are the non-resonant terms, which 
have no analogy in the non-relativistic model. The third term can be reduced (jp 2 — ► — p 2 ) 
to a similar term in the non-relativistic model stemming from the A-hole excitation, but 
there exist some differences mainly caused by the different masses of nucleons and deltas. 
The situation might be understood in view that the problem of describing a spin 3/2 
particle in relativistic quantum field theory remains unsolved. Fortunately, the difference 
between the non-relativistic limit of the relativistic model and the standard non-relativistic 
model is quantitatively insubstantial. 



Fig. 12 



Fig. 12(a) displays the pion dispersion relation (the pion branch) at normal density. 
The solid line denotes the free pion dispersion relation. The dotted line is calculated with 
Eqs. (E5) and (E6), i.e., the non-relativistic limit of the relativistic model but with the 
relativistic kinetics. The dashed line is computed by taking p 2 — > — p 2 in Eqs. (E5) and 
(E6). One can clearly see that the relativistic effect (mainly from the kinematic origin) 
makes the pion dispersion relation harder at low momenta and softer at high momenta. 
Furthermore, the relativistic effect at low momenta mainly comes from the non-resonant 
terms, i.e., the first and second terms on the r.h.s. of Eq. (E6). If one switches off these 
two terms, at low momenta the results (the dash-dotted line) approach to the dashed line 
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while at high momenta to the dotted line. Fig. 12(b) depicts the pion dispersion relation 
(both the pion and the A-hole branches) at different densities calculated with Eqs. (E5) 
and (E6) and p 2 —>■ — p 2 . The short-range correlation effect is included through Eq. 
(125) (p 2 — > — p 2 ) with g' = 0.6. As one can see from the figure, the obtained pion 
dispersion relation is nearly the same as that of the non-relativistic model []T6"|,|2"T[| except 



the pion branch is a little harder at high momenta. In this case one may conclude that 
the difference between the results of Fig. 6 and that of the non-relativistic model mainly 
stems from the relativistic kinetics. 

APPENDIX F 

In this appendix we derive the conserved current and energy-momentum tensor. 
(1) Current 

Make a four-dimension integration of momentum on the both side of Eq. (62), the 



right-hand side (the collision term) goes to zero [77] and we have 



^iten^p(x, P )ag} /ir ^ r) = o. (Fi) 



It is straightforward to find the current conservation 

= (F2) 

with 

/ w (x,p,r) 



J M (x) 



(F3) 



u*(p) 

We note that each /^(x, p, r) is in principle accompanied by a 5 function 5(p — w*(p)) 

for on-shell pion. 

(2) Energy-momentum tensor 

Multiply Pn on both side of Eq. (62) and make a four-dimension integration of mo- 
mentum we arrive at 
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(F4) 



Our strategy is to extract the d v x out of the whole equation. For the first term it is 
straightforward. The second term can be rewritten as 



2lu*(p) 

With the help of Eqs. (73) and (74), it becomes 



3: 



The third and fourth terms can be written as 



Jd 4 pp,dA\ ( d * ReU ^(x,P)) 



A(x,p,r) 

u*(p) 



J d 4 pp,d: l -(dlReYlT~ p {x,p)) 



A(x,p,r) 

u*(p) 



= d A p\ (d^ReU^p)) - er x J d 4 p\p, (d^ReU r - p (x,p)) 



-d:Jd 4 p^p,(d^ReUr o - p (x,p 



A(x,p,r) 



(F5) 



The second term of the above equation turns out to be 

= d 4 PPx (d x p Reu ro ~ p (x,p)) % fvi *;?' T) + ^M^(x )P )^r Wx ' p,r) 



u)*{p) 

/tt(x,P,t ) , 1 

cu*(p) 



1 (' 

J d 4 pp x (d^ReUj-- p (x,p)) d[ 

d:\g„J d 4 PPx (diReu r - p (x, P )) u *;* p) T) . 



J d 4 pp x (d^ReU loop (x,p 



LJ*(p) 

/^(x,p,r) 

u*{p) 

(F6) 



In the first and second equality of Eq. (F6) we have used the fact that the terms with the 
double derivative acting on the same quantity can be neglected in the gradient expansion. 
At the end we have the energy-momentum conservation 



d» x T^{x) = 



(F7) 
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with 

- g^\j d 4 pReU^(x,p)^^^- 

- \ J d 4 PP, (d^ReU T - p (x,p)) Ui *;V p) T) 

+ J d 4 PPx (d x p ReU r - p (x,p)) fA ^l T) . (F8) 
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TABLES 



TABLE I: Symbols and notation used for the baryon-baryon-meson vertex, is the 
transformed four-momentum. 



A 


m A 


&NN 


Saa 


7a 




T A 


$ A (x) 


^A 


A 


a 


m CT 


&NN 


Saa 


1 


1 


1 


<j(x) 


1 


1 


UJ 


m w 




Saa 




1 


1 






1 


7T 




IT 

&NN 


7T 

Saa 


f75 


T 


T 


7T(x) 


1 




P 


m p 


2&NN 


2 Saa 


7m 


T 


T 


p^(x) 


_ g/f 


<5« 



TABLE II: Symbols and notation used for the interaction vertex involving only mesons. 
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07T7T 




r A 


$ A (x) 
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2"So"7T 


1 




cr(x) 
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gpn 






Pk( x ) 
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CAPTIONS 



Fig.l Contour along the axis for an evaluation of the operator expectation value. In 
practice, to is shifted to — oo and t max to +00. 

Fig. 2 Feynman diagrams contribute to the Hartree term of pion self-energy. A wavy 
line denotes the exchanged virtual meson, a solid line, a doubled line and a dashed 
line represents the nucleon, delta and pion, respectively. 

Fig . 3 Feynman diagrams contribute to the Fock term (one baryon-loop) of pion self- 
energy. Different lines denote the different particles as described in the caption of 
Fig. 2. 

Fig . 4 Feynman diagrams contribute to the Born term of pion self-energy. Different lines 
denote the different particles as described in the caption of Fig. 2. The imaginary 
part of Fig. 4(a) contributes to the 71N — > ttN elastic cross section, and Fig. 4(b) 
to the it A — > 7rA, Fig. 4(c), (d) to the irir — > -kit elastic cross section, respectively. 

Fig. 5 The gap between the effective masses of particles and anti-particles as a function 
of density. The universal coupling strengths for the nucleons and deltas are assumed. 

Fig. 6 The pion dispersion relation in symmetric nuclear matter at saturation density. 
The solid line shows the free dispersion relation. The dotted curves represent the 
in-medium dispersion relation for different branches as indicated in the figure. The 
upper and lower hatched areas indicate regions of non- vanishing imaginary parts of 
the A-hole and nucleon-hole polarizations, respectively. 

Fig. 7 (a) The pion dispersion relation at different densities, (b) The real part of pion 
self-energy at different densities relating to the pion dispersion relation shown in 
(a), (c) The contributions of different excitation modes to the real part of pion 
self-energy. The calculations are performed at normal density. 
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Fig. 8 The pion dispersion relation (dotted curves) at different densities. The short- 
range correlations are taken into account in a non-relativistic way with Migdal 
parameter g'=0.6. The solid lines represent the free dispersion relation. 

Fig. 9 Free cross sections for reactions ir + p — > A ++ and ir~p — > A . The dots are 
the experimental data from Ref. [E3]. The solid curve is our results calculated 



with the mass-dependent coupling strength while the dashed curve with the mass- 
independent coupling strength. 

Fig. 10 The momentum-dependent A-decay width in the free space. The solid line and 
long dash-dotted line are our results computed with MDC and MIC, respectively. 
Dotted pl |, dashed |85| and dash-dotted line represent several phenomenological 



parameterizations commonly used in the transport models. 

Fig. 11 The in- medium A-formation cross section and A-decay width. In (a) and (b) the 
free pion mass is used in the calculations. In (c) and (d) the effective pion mass is 
taken into account via the dispersion relation of Fig. 7(a). 

Fig. 12 (a)The pion dispersion relation (the pion branch) at normal density. Different 
lines correspond to the different situations as explained in text. (b)The pion disper- 
sion relation (the pion and A-hole branch) at different densities. The calculations 
are performed with Eqs. (E5) and (E6) and p 2 ^ — > — p 2 . The short-range correlation 
effect is taken into account by means of the Migdal parameter of g' = 0.6. 
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Fig. 2 
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Fig. 5 
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Fig. 8 
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Fig. 10 
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Fig. 11 




65 




66 



